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PREFACE 


This  report  is  part  of  Rand's  continuing  study  of  the  effects  of 
abnormally  ionized  atmospheres,  such  as  those  produced  by  nuclear  bursts, 
on  the  propagation  of  radio  waves  and  on  military  radio  communication 
systems,  present  and  future.  As  a  tool  useful  for  providing  answers 
to  many  problems  arising  in  such  studies,  we  have  developed  a  general 
computer  program  which  could  be  applied  over  a  wide  range  of  frequen- 
cies  (typically  10  to  10  Hz)  for  any  specified  vertical  ionization 
structure  of  the  atmosphere.  Both  electrons  and  positive  and  negative 
ions  are  included  in  the  program.  For  specified  frequencies,  angle  of 
wave  incidence  on  the  ionosphere,  and  ionization  profiles  the  program 
provides  the  real  and  imaginary  components  of  the  complex  refractive 
index  at  all  heights,  the  local  wavelength  of  the  signal,  the  differ¬ 
ential  absorption  at  all  heights,  and  the  reflection  and  transmission 
parameters  for  the  nonpenetrating  and  penetrating  signal  components, 
respectively.  The  program  provides  results  for  either  quasi-longitu¬ 
dinal  or  quasi-transverse  conditions  of  transmission. 

The  program  is  developed  from  a  ray  theory  approach  and  as  such 
has  certain  restrictions  and  accuracy  limitations,  in  particular  at 
the  lower  frequencies  used.  Despite  such  limitations  the  program  is 
a  very  versatile  tool  for  obtaining  results  which  can  be  applied  to 
practical  assessment  of  many  problems  in  radio  transmission,  both  for 
earth-to-earth  and  earth-to-space  radio  links,  and  for  ambient  or  arbi¬ 
trarily  disturbed  ionospheres. 

Selected  results  obtained  from  the  program  are  presented  in  this 
and  a  companion  report,  R-558-PR,  Transmission  of  Electromagnetic  Waves 
Through  Normal  and  Disturbed  Atmospheres . 

Henry  G.  Booker  is  Chairman  of  the  Department  of  Applied  Physics 
and  Information  Science  at  the  University  of  California,  San  Diego  and 
is  a  consultant  to  The  Rand  Corporation. 
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SUMMARY 


This  report  presents  the  technical  basis  used  in  the  development 

of  a  computer  program  for  assessing  the  effects  of  ambient  and  disturbed 

daytime  and  nighttime  ionospheres  on  the  reflection  and  transmission 

10  -3 

of  radio  waves  in  the  10  to  10  Hz  frequency  range.  The  ionospheric 
models  and  electromagnetic  theory  used  are  described,  certain  limita¬ 
tions  are  discussed,  and  the  parameters  and  equations  used  are  summa¬ 
rized.  A  number  of  results  obtained  from  the  program  are  presented  in 
graphical  form.  From  these,  the  existence  of  height-dependent  pass 
bands,  stop  bands,  and  conduction  bands  in  the  ionosphere  is  evident. 

Also  shown  are  reflection  levels  and  two-way  reflection  losses.  These 
results  provide  useful  understanding  of  ionospheric  wave  propagation, 
and  serve  as  background  material.  Results  with  direct  engineering  ap¬ 
plication  to  earth-satellite  communication  links  have  also  been  obtained, 
and  are  given  in  a  companion  report. 
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I .  INTRODUCTION 


In  this  report  a  technical  description  of  a  computer  program  we 
have  developed  is  given,  along  with  some  diagrams  which  present  a  use¬ 
ful  overall  impression  of  the  effects  of  the  ionosphere  on  both  reflec¬ 
tion  and  transmission  over  the  entire  range  of  frequencies  from  micro- 

10  -3 

wave  to  hydromagnetic  (10  to  10  Hz).  Ambient  day  and  night  condi¬ 
tions  are  treated,  as  well  as  disturbed  conditions  corresponding  either 
to  a  mild  nuclear  environment  or  a  severe  polar  cap  absorption  (PCA) 
event.  In  addition  to  the  material  included  in  this  report,  which  is 
useful  for  purposes  of  general  interpretation  and  understanding  of 
propagation  characteristics,  we  have  obtained  results  with  direct  engi¬ 
neering  applicability  to  satellite-earth  communication  links.  These 
are  presented  in  a  companion  report ^  in  the  form  of  coefficients  of 
transmission  through  the  ionosphere  and  height  profiles  of  attenuation. 

Section  II  describes  the  assumptions  made  and  the  models  used, 
and  Section  III  gives  the  numerical  results  and  related  discussion. 

The  Appendix  details  the  wave-propagation  equations  and  the  parameters 
used  in  the  computer  program  from  which  the  results  of  the  report  were 
derived.  The  wave  equations  are  those  of  the  magneto-ionic  theory, 
taking  both  electrons  and  ions  into  account. 
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II.  MODEL  AND  ASSUMPTIONS 


Models  of  the  ionosphere  are  required  that  specify  as  functions 
of  height  the  electronic  and  ionic  collisional  frequencies,  the  elec¬ 
tron  density,  and  the  densities  of  positive  and  negative  ions.  Dif¬ 
ferent  models  are  needed  for  night  and  day.  Also  of  interest  are  con¬ 
ditions  in  which  the  D  region  of  the  ionosphere  is  disturbed  by  natural 
or  man-made  phenomena.  The  four  profiles  of  electron  and  negative  ion 
density  used  in  this  study  are  shown  in  Fig.  1  and  are  identified  as 
ambient  night,  ambient  day,  disturbed  night,  and  disturbed  day. 

The  ambient  profiles  of  Fig.  1  are  typical  of  those  in  current 
use.  Positive  ion  profiles  were  chosen  to  maintain  charge  neutrality. 
The  collisional  frequencies  used  are  based  on  Ref.  2.  The  principal 
uncertainty  in  the  profiles  involves  the  distributions  of  electron 

g 

density  below  a  value  of  about  10  electrons  per  cubic  meter,  where 
measurements  are  the  most  uncertain.  The  profiles  for  the  disturbed 
day  and  night  cases  are  based  on  electron  production  rates  due  to  wide¬ 
spread  high-altitude  fission  debris  as  described  in  Ref.  3,  using  a 

* 

value  of  5000  for  ,  and  nominal  reaction  rates  such  as  given  in 
Ref.  4.  Use  of  production  and  reaction  rates  as  recently  given  in 
Ref.  5  would  be  expected  to  lead  to  better  prediction  of  the  electron 
and  ion  distributions;  however,  the  differences  are  small  and  not  sig¬ 
nificant  for  the  purposes  of  this  report.  Although  the  disturbed  ion¬ 
ization  profiles  are  representative  of  a  mild  nuclear  environment  from 
high-altitude  fission  debris,  the  values  of  electron  and  ion  density 
in  the  important  height  region  from  about  40  to  80  km  are  approximately 
those  expected  for  a  severe  polar-cap  disturbance,  and  hence  the  trans¬ 
mission  and  reflection  results  can  be  taken  to  be  appropriate  for  such 
a  disturbance.  The  ambient  and  disturbed  cases  thus  represent  the  ap¬ 
proximate  lower  and  upper  limits  of  naturally  occurring  ionospheric 
structures . 

The  major  approximation  is  the  use  of  WKB  methods  to  compute|the 
transmission  and  reflection  losses.  This,  of  course,  neglects  the 

*P  c  lr)13  fission  yield  in  megatons 

M  (A  u  •  AA  •  i  \2  „1.2 

(debris  radius  m  km)  t 

sec 


Height  (kilometers) 


Charged-particle  density  (per  cubic  meter) 

Fig.  1  —  Ionospheric  models  for  the  ambient  night,  ambient  day,  disturbed  night, 

and  disturbed  day  ionospheres 
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effects  of  gradient  reflection  and  is  valid  provided  that  the  wave¬ 
length  in  the  ionosphere  is  smaller  than  the  characteristic  distance 
over  which  the  ionospheric  refractive  index  changes  substantially — 
typically  a  few  kilometers.  In  Fig.  2  we  show  the  ratio  of  the  wave¬ 
length  in  the  ionosphere  to  the  free-space  wavelength  as  a  function 
of  height  for  a  range  of  frequencies  and  ambient  day  and  night  condi¬ 
tions.  The  ratios  were  computed  from  the  equations  given  in  the  Ap¬ 
pendix.  Note  that  the  wavelength  is  shortened  considerably  in  the 
ionosphere,  especially  at  the  lower  frequencies.  By  multiplying  the 
curves  of  Fig.  2  by  the  appropriate  free-space  wavelength  to  get  the 
height-dependent  absolute  wavelength,  it  can  be  seen  that  the  WKB  so¬ 
lutions  should  be  quite  good  for  frequencies  above  about  1  kHz.  Be¬ 
low  about  100  Hz  the  WKB  validity  condition  is  seriously  violated. 

Thus,  at  frequencies  of  the  order  of  100  Hz  or  less,  gradient  reflec¬ 
tion  is  not  negligible  and  our  results  for  attenuation  correspond  to 
lower  bounds  on  the  transmission  losses.  Similarly,  resonances  in  the 
transmission  coefficients  and  Brewster's  angle  effects  in  reflection 
coefficients  cannot  be  correctly  predicted  by  WKB  methods  at  these  low 
frequencies.  We  have  spot  checked  some  of  the  results  at  frequencies 
of  the  order  of  1  Hz  by  using  a  full-wave  numerical  method  developed 
by  Rand  for  the  study  of  the  ionospheric  transmission  of  geomagnetic 
micropulsations . In  cases  where  meaningful  comparisons  could  be 
made,  the  agreement  was  within  an  order  of  magnitude.  This  indicates 
that  even  at  micropulsation  frequencies  the  WKB  methods  yield  useful, 
albeit  rough,  estimates.  Our  conclusion  that  the  WKB  methods  should 
be  very  reliable  at  frequencies  above  1  kHz  is  confirmed  by  the  work 
of  Jones, who  compared  full-wave  and  WKB  calculations  of  ionospheric 
reflection  of  waves  in  the  1  kHz  to  10  MHz  band.  In  no  case  were  the 
WKB  calculations  in  error  by  more  than  a  few  decibels. 

The  calculations  have  been  performed  for  the  two  cases  of  waves 
having  phase  velocities  in  the  vertical  direction  and  at  an  angle,  i, 
of  80  deg  with  respect  to  the  vertical  at  the  base  of  the  ionosphere. 
The  geomagnetic  field  was  assumed  vertical.  Thus,  strictly  speaking, 
the  results  for  an  incidence  angle,  i,  of  zero  deg  correspond  to  verti¬ 
cal  propagation  at  a  magnetic  pole,  i.e.,  vertical  longitudinal 
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propagation.  Booker  and  Dycev  have  investigated  the  range  of  validity 
of  the  quasi-longitudinal  (QL)  approximation  for  hydromagnetic  waves 
and  radio  waves  in  the  ionosphere  and  magnetosphere  in  the  absence  of 
collisions.  We  have  here  re-examined  the  validity  range  after  includ¬ 
ing  collisions,  which  increase  the  accuracy  of  the  QL  approximation, 
and  for  vertical  propagation  the  conclusions  are  these: 

1.  At  altitudes  above  about  150  km,  where  electron-neutral  and 
ion-neutral  collisions  tend  to  be  relatively  unimportant,  the 
QL  approximation  is  valid  for  geomagnetic  dip  angles  larger 
than  40  deg  if  the  wave  frequency  exceeds  about  10  Hz.  For 
wave  frequencies  less  than  about  1  Hz,  the  QL  approximation 
applies  only  near  the  magnetic  pole. 

2.  At  an  altitude  of  120  km,  where  ion-neutral  collisions  are 
important,  the  QL  approximation  is  valid  for  dip  angles  larger 
than  40  deg  for  essentially  all  frequencies  considered  here. 

3.  At  altitudes  below  about  70  km,  where  electron-neutral  colli¬ 
sions  dominate,  the  QL  approximation  is  valid  everywhere. 

Note  that  at  the  altitudes  where  the  most  important  absorption 
occurs  (below  150  km),  the  QL  approximation  applies  down  to  mid-lati¬ 
tudes  or  lower.  A  detailed  application  of  the  QL  approximation  would 
involve  recalculating  the  results  to  be  shown  using  the  equations  of 
the  Appendix,  but  at  each  latitude  substituting  the  vertical  component 
of  the  geomagnetic  field  for  the  total  field  which  now  appears  in  the 
formulas.  This  has  been  done  in  some  cases  as  a  check  and,  for  lati¬ 
tudes  above  about  40  deg,  the  losses  were  only  somewhat  larger  than 
those  obtained  from  the  vertical  polar  calculations,  the  results  of 
which  are  presented  in  Section  III. 

At  frequencies  less  than  the  electronic  gyro  frequency,  the  direc¬ 
tion  of  group  propagation  (the  ray)  can  differ  drastically  from  the 
direction  of  phase  propagation  (the  wave  normal)  even  within  the  domain 
of  validity  of  the  QL  approximation.  (See  Ref.  8  for  a  detailed  analy¬ 
sis  of  this  phenomenon  in  the  absence  of  collisional  effects.)  At  VLF 
frequencies  and  lower,  the  ionospheric  refractive  index  is  large,  and 
a  wave  traveling  downward  from,  say,  a  satellite  will  have  its  wave 
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nonnal  bent  strongly  away  from  the  vertical  upon  emergence  from  the 
ionosphere.  (Actually  this  occurs  over  a  height  range  in  the  lower 
ionosphere.)  Conversely,  a  low-frequency  wave  incident  obliquely  upon 
the  ionosphere  from  below  will  have  its  wave  normal  bent  to  nearly  the 
vertical  direction  within  a  few  ionospheric  conductivity-scale  heights. 
In  order  to  emerge,  rather  than  be  internally  reflected,  a  downward 
propagating  low-frequency  wave  must  have  its  normal  very  nearly  verti¬ 
cal  while  in  all  but  the  lowest  regions  of  the  ionosphere.  Thus  our 
results  given  here  and  in  Ref.  1,  which  correspond  to  vertical  or 
nearly  vertical  phase  propagation  in  the  ionosphere ,  apply  to  realistic 
satellite-ground  phase  paths.  This  is  true  even  though  the  group  path 
tends  to  be  more  or  less  along  the  geomagnetic  field  in  the  absence  of 
collisions'1  and,  for  temperate  latitudes,  can  be  distinctly  non¬ 
vertical.  The  relation  between  ray  and  phase  propagation  direction  is 
not  well  understood  when  both  collisional  and  geomagnetic  effects  are 
of  consequence.  In  any  event,  however,  the  transmission  loss  is  ob¬ 
tained  by  a  vertical  or  nearly  vertical  integration  over  all  but  the 
lowest  ionospheric  heights  and  our  results  should  in  most  cases  provide 
reasonable  estimates. 

As  indicated  earlier,  results  also  will  be  given  for  glancing  in¬ 
cidence  (or  emergence),  i  =  80  deg.  However,  the  longitudinal  assump¬ 
tion  is  retained  for  the  glancing-incidence  calculations.  Thus  the 
curves  to  be  shown  for  glancing  incidence  give  only  an  indication  of 
where  angle  of  incidence  matters  and  of  roughly  how  much.  Recalling 
that  a  typical  phase  path  from  a  VLF/ELF  satellite  transmitter  will 
be  nearly  vertical  until  the  lower  ionosphere  is  reached,  whence  the 
propagation  vector  is  bent  to  nearly  horizontal,  it  is  evident  that  for 
certain  satellite  -  ground-station  geometries  a  portion  of  the  phase 
propagation  path  could  be  nearly  perpendicular  to  the  geomagnetic  field, 
even  at  high  latitudes.  For  such  paths,  a  full  analysis  involving 
neither  QL  nor  QT  approximations  is  needed.  Such  aspects  are  beyond 
the  scope  of  the  treatment  given  in  this  report.  It  is  worth  mention¬ 
ing,  however,  that  we  have  performed  computations  for  the  case  where 
phase  propagation  is  transverse  to  the  geomagnetic  field.  We  will  not 
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present  detailed  graphs  of  these  results  here  but  will  simply  point 
out  that  the. attenuation  for  transverse  propagation  was  found  to  be 
larger  than  for  longitudinal  propagation  if  all  other  variables  are 
held  constant.  This  indicates  that  if  the  QL  approximation  is  vio¬ 
lated  over  part  of  a  phase  path,  the  attenuation  and  transmission 
losses  will  be  larger  than  given  here  or  in  Ref.  1. 

Finally,  note  that  the  effects  of  earth  curvature  are  neglected. 


-9- 


III.  RESULTS  AND  DISCUSSION 


-3 

Figures  3  to  10  show  frequency  plotted  horizontally  from  10  to 
10  Hz.  They  show  ionospheric  height  plotted  vertically,  with  heights 
above  180  km  omitted,  since  it  is  reasonably  obvious  what  happens  beyond 
that  point.  Significant  frequencies  in  the  figures  are  the  following: 

1.  The  superior  penetration  frequencies,  or  "critical  frequen¬ 
cies,"  f  and  fv„  of  the  F  region  of  the  ionosphere  for  the 
0  and  X  waves  (that  is,  for  left-handed  and  right-handed  po¬ 
larization  about  the  earth's  magnetic  field,  respectively). 

2.  The  electronic  gyrofrequency  f  . 

3.  The  ionic  gyrofrequency  f  separated  into  an  atomic  value 
f.,  above  150  km  and  a  molecular  value  f„  below  160  km,  with 
an  overlap  of  atoms  and  molecules  between  these  two  heights . 

(In  the  model  positive  ions  were  considered  molecular  below 
150  km  and  atomic  above  160  km.  The  negative  ion  density 
above  90  km  was  assumed  zero.) 

Figures  3  to  6  refer  to  night  models  and  Figs.  7  to  10  to  day 
models.  Alternate  figures  refer  to  the  0  and  X  waves.  Figures  3  and 
4  refer  to  ambient  (undisturbed)  conditions  at  night),  while  Figs.  7 
and  8  refer  to  ambient  conditions  in  the  daytime. 

In  each  diagram  is  shown  the  curve  for  which 

ve  =  Ke  1  “I  (1) 

where  is  the  electronic  collisional  frequency,  wMe  is  the  angular 
electronic  gyrofrequency,  id  is  the  angular  wave  frequency,  and  the 
upper  and  lower  signs  refer  to  the  0  and  X  waves,  respectively.  At 
heights  above  this  curve  the  effect  of  electronic  collisions  is  rela¬ 
tively  secondary.  Below  the  level  where  ■=  (between  70  and  80  km) 
the  effect  of  the  earth's  magnetic  field  is  unimportant.  At  an  alti¬ 
tude  between  110  and  120  km  the  level  for  which  v.  =  ul,.  can  also  be 

l  Mi 

seen,  where  is  the  ionic  collisional  frequency  and  is  the  ionic 
gyrofrequency.  Above  this  level  both  ionic  and  electronic  collisions 
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Fig.  3 — Quasi-longitudinal  propagation 


of  the  0  wave  for  ambient 


(undisturbed)  night  conditions 
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Fig.  5 — Quasi-longitudinal  propagation  of  the  0  wave  for  disturbed  night  conditions 


Reflection  loss  (decibels) 


Fig.  6 — Quasi-longitudinal  propagation  of  the  X  wave  for  disturbed  night  conditions 
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Fig.  7— Quasi-longitudinal  propagation  of  the  0  wave  for  ambient 


(undisturbed)  day  conditions 
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Fig.  8 — Quasi-longitudinal  propagation  of  the  X  wave  for  ambient  (undisturbed)  day  conditions 
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Fig.  9  Quasi-longitudinal  propagation  of  the  0  wave  for  disturbed  day  conditions 


Height  (kilometers) 


Fig.  10 — Quasi-longitudinal  propagation  of  the  X  wave  for  disturbed  day  conditions 


Reflection  loss  (decibels) 
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are  secondary,  although  they  may  still  be  responsible  for  some  inter¬ 
esting  attenuation. 

Above  the  level  where  we  see  in  Figs.  3  to  10  the  famil¬ 

iar  magneto-ionic  behavior  for  quasi-longitudinal  propagation  in  the 
absence  of  collisions.  This  involves  a  pass  band  below  the  ionic  gyro- 
frequency  for  both  the  0  and  X  waves  (the  hydromagnetic  pass  band), 
and,  between  the  two  gyrofrequencies,  a  stop  band  for  the  0  wave  and 
a  pass  band  for  the  X  wave  (the  whistler  pass  band) .  Above  the  usual 
critical  frequency  there  is  the  electromagnetic  pass  band.  In  Figs.  3 
to  10,  pass  bands  are  unstriped  and  stop  bands  are  striped  horizontally. 

It  will  be  noticed  in  Figs.  3  to  10  that  the  behavior  above  the 
level  vi  =  between  the  two  gyrofrequencies  extends  below  the  ionic 
gyrofrequency  in  the  altitude  range  between  vi  =  and  . 

This  is  true  for  both  the  0  and  X  waves.  In  Fig.  3  the  whistler  stop 
band  for  the  0  wave  extends  down  through  the  hydromagnetic  frequency 
band  in  the  height  range  between  v  =  or,,  and  v  =  ;  whether  or  not 

this  stop  band  below  the  ionic  gyrofrequency  is  thick  enough  in  rela¬ 
tion  to  the  wavelength  in  the  medium  to  be  effective  as  a  stop  band 
will  be  investigated  later  in  this  report.  In  Fig.  4  the  whistler 
pass  band  for  the  X  wave  extends  down  through  the  hydromagnetic  fre¬ 
quency  band  in  the  height  range  between  v.  =  uw.  and  v  =  to  .  thereby 

l  Mi  e  Me 

converting  the  entire  region  for  which  —  to  into  a  single  pass 

band  for  the  X  wave.  Whistler  behavior  extends  below  the  ionic  gyro¬ 
frequency  in  the  height  interval  between  v.  =  to...  and  v  =  to,,  because, 
in  these  circumstances,  collisions  inhibit  motion  of  the  ions  while 
permitting  motion  of  the  electrons.  This  is  what  happens  for  a  dif¬ 
ferent  reason  under  normal  whistler  conditions,  where  it  is  the  com¬ 
bined  effect  of  inertia  and  the  earth's  magnetic  field  that  inhibits 
ion  and  permits  electron  motion. 

Besides  pass  bands  and  stop  bands,  there  is  shown  at  sufficiently 
low  frequencies  in  each  of  Figs.  3  to  10  a  conduction  band  indicated 
by  vertical  striping.  This  is  the  region  where  the  square  of  the  re¬ 
fractive  index  is  dominated  by  its  imaginary  part,  so  that  the  refrac¬ 
tive  index  has  properties  similar  to  those  encountered  in  a  metal.  The 
thickness  of  the  conduction  band  varies  with  frequency,  and  this 
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thickness  has  to  be  compared  with  the  "skin  depth"  before  we  know 
whether  the  conduction  band  is  mainly  transparent  or  mainly  reflect¬ 
ing  at  any  particular  frequency. 

The  pass,  stop,  and  conduction  bands  in  Figs.  3  to  10  for  verti¬ 
cal  propagation  may  be  defined  as  the  regions  where  the  square  of  the 
refractive  index  is  mainly  real  and  positive  (pass  band),  mainly  real 
and  negative  (stop  band),  and  mainly  imaginary  and  negative  (conduc¬ 
tion  band).  Following  Crain  and  Booker, ^  we  may  take  the  refractive 
index  n  as  given  by 


n2 3  =  1  -  A  -  jB  (2) 

* 

where  A  and  B  are  real,  and  define  boundaries  for  angle  of  incidence 
i  as  follows: 

1.  The  boundary  between  a  pass  band  and  a  stop  band  may  be 
taken  as 

2 

A  =  cos  i  (3) 

2.  The  boundary  between  a  conduction  band  and  a  pass  band  may 
be  taken  as 

B  =  cos2i  (4) 

3.  The  boundary  between  a  stop  band  and  a  conduction  band  may 
be  taken  as 


A  =  B  (5) 

(9) 

As  shown  by  Field  and  Engel  a  coefficient  somewhat  greater  than 
unity  on  the  right-hand  side  of  Eq.  (4)  is  desirable.  Such  a  coeffi¬ 
cient  was  in  fact  used  in  the  calculations;  the  curves  showing  the 

* 

See  the  Appendix  for  equations  which  express  A  and  B  in  terms 
of  ionospheric  properties. 


-20- 


low-frequency  boundary  of  the  electromagnetic  pass  band  then  require 
slight  smoothing  where  the  phenomenon  on  the  other  side  of  the  boundary 
changes  from  conduction-band  behavior  to  stop-band  behavior. 

The  low-frequency  boundary  of  the  electromagnetic  pass  band  is 
given  by  Eq.  (3)  when  vg  <  [  ±  oj  [  and  by  Eq .  (4)  when  >  |  ±  to  |  ; 

therefore  it  depends  on  angle  of  incidence.  The  striped  regions  shown 
in  Figs.  3  to  10  refer  to  vertical  incidence  (i  =  0  deg).  The  extent 
to  which  these  should  be  extended  towards  higher  frequencies  for  glanc¬ 
ing  incidence  is  indicated  by  the  almost  parallel  curve  marked  i  =  80  deg. 

At  high  frequencies  the  edge  of  the  electromagnetic  pass  band 
shown  in  Figs.  3  to  10  is  usually  regarded  as  the  height  of  reflection. 

In  fact,  of  course,  reflection  takes  place  from  a  stratum  extending 
above  this  height  to  a  level  conveniently  defined  as  the  level  where 
the  storage  field  has  dropped  to  about  1  neper  below  the  value  at  the 
height  of  reflection.  Using  the  WKB  approximation,  this  level  for 
1  neper  penetration  has  been  calculated  and  plotted  in  Figs.  3  to  10. 

In  Fig.  3,  the  reflecting  stratum  acquires  more  than  "pencil  thick¬ 
ness"  below  about  10“*  Hz  for  vertical  incidence  and  10^  Hz  for  glancing 
incidence.  The  curve  for  1  neper  penetration  intersects  the  top  of  the 
stop  band  just  below  the  ionic  gyrofrequency .  At  all  lower  frequencies 
the  combined  stop  and  conduction  band  for  the  0  wave  has  a  thickness 
less  than  the  "skin  depth"  as  understood  for  a  metal.  Below  about  10  Hz 
in  Fig.  3,  therefore,  the  conduction  and  stop  bands  are  in  fact  trans¬ 
parent  . 

The  frequency  at  which  the  curve  for  1  neper  penetration  inter¬ 
sects  the  top  of  the  stop  band  for  the  0  wave  (Figs.  3,  5,  7,  9)  may 
be  defined  as  the  hydromagnetic,  or  inferior,  penetration  frequency  of 
the  ionosphere  for  the  0  wave.  A  similar  inferior  penetration  frequency 
of  the  ionosphere  can  be  defined  for  the  X  wave  by  the  intersection  of 
the  curve  for  1  neper  penetration  with  the  top  of  the  conduction  band 
(Figs.  6  and  10).  In  the  case  of  Fig.  4  the  attenuation  in  the  conduc¬ 
tion  band  at  vertical  incidence  is  less  than  1  neper,  so  that  the  con¬ 
duction  band  is  to  be  classified  as  transparent;  it  acts  like  a  metal 
film  whose  thickness  is  small  compared  with  the  skin  depth.  The  inter¬ 
pretation  of  Fig.  8  requires  care  and  is  discussed  later. 
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Between  the  inferior  penetration  frequency  of  the  ionosphere  just 

discussed  and  the  usual,  or  superior,  penetration  frequency  known  as 

k 

the  critical  frequency,  the  ionosphere  reflects  .  Nevertheless  the 
attenuation  of  the  reflected  wave  below  the  level  of  reflection  may 
be  quite  substantial.  This  attenuation  is  therefore  shown  in  Figs.  3 
to  10,  using  the  right-hand  vertical  scale. 

The  attenuation  below  the  level  of  reflection  in  Figs.  3  to  10 
has  been  calculated  by  the  methods  developed  by  Booker  and  Crain^^ 
and  by  Booker,  Fejer,  and  Lee.  In  view  of  the  criticism  by  von 

Roos^^  of  the  theorem  given  by  Booker,  Fejer,  and  Lee,  ^ ^ ^  it  should 
perhaps  be  emphasized  that  the  theorem  does  in  fact  apply  with  adequate 
accuracy  at  all  levels  in  Figs.  3  to  10  below  the  curves  corresponding 
to  Eq.  (1).  Von  Roos  is  quite  correct  in  stating  that  the  theorem  is 
unsatisfactory  when  the  effect  of  the  earth's  magnetic  field  is  vital. 
Thus,  it  is  unsatisfactory  for  calculating  the  attenuation  of  the  re¬ 
flected  0  wave  in  the  VHF  band  under  ambient  night  conditions  (Fig.  3); 
however,  the  attenuation  is  then  relatively  small. 

Figures  3  to  10  show  superior  penetration  frequencies  (critical 
frequencies)  for  the  ionosphere  that  depend  on  maximum  electron  density 
and  angle  of  incidence  (based  on  a  plane  earth  and  ionosphere)  in  the 
usual  way.  On  the  other  hand,  they  also  indicate  inferior  penetration 
frequencies  for  the  ionosphere  arising  from  thinness  of  the  reflecting 
stratum  as  listed  in  Table  1. 

Of  course,  the  inferior  penetration  frequency  of  the  ionosphere  is 
not  as  precise  a  concept  as  the  superior  penetration  frequency.  In  the 
neighborhood  of  the  inferior  penetration  frequency  there  is  an  interval 
of  about  an  octave  over  which  significant  partial  reflection  and  pene¬ 
tration  is  taking  place.  However,  there  is  a  sense  in  which  the  inferior 
penetration  frequency  tends  to  be  more  precise  than  the  superior  pene¬ 
tration  frequency.  The  inferior  penetration  frequency  shows  little 
dependence  on  angle  of  incidence  because,  at  the  inferior  penetration 

k 

This  reflection  occurs  within  the  context  of  the  WKB  treatment 
used  here.  As  discussed  earlier,  a  full-wave  treatment  would  indicate 
significant  gradient  reflection  at  frequencies  of  100  Hz  or  less. 
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Table  1 

INFERIOR  PENETRATION  FREQUENCY  OF  THE  IONOSPHERE  IN  HERTZ 
(QUASI-LONGITUDINAL  PROPAGATION) 


0  Wave 

X  Wave 

Condition 

O 

O 

II 

*H 

e 

O 

00 

II 

H* 

II 

O 

o 

i  =  80° 

Ambient  night 

fMi 

fMi 

fMe 

fMe 

Disturbed  night 

0.6 

0.6 

200 

300 

Ambient  day 

1.6 

1.6 

io4 

104 

Disturbed  day 

0.4 

0.4 

25 

25 

frequency,  the  refractive  index  in  the  ionosphere  is  usually  high,  so 
that  waves  incident  at  almost  any  angle  in  free  space  are  refracted  to 
nearly  vertical  propagation  in  the  ionosphere.  The  inferior  penetra¬ 
tion  frequencies  shown  in  Table  1  are  reasonably  close  to  the  maximum 
frequencies  of  the  lower  ionospheric  transmission  bands  defined  on  the 
basis  of  a  different  criterion  in  Ref.  1. 

Table  1  shows  that  under  ambient  night  conditions,  the  inferior 
penetration  frequency  is  given  essentially  by  the  ionic  gyrofrequency 
for  the  0  wave  and  by  the  electronic  gyrofrequency  for  the  X  wave. 
Collisions  are  relatively  unimportant  at  any  level  where  appreciable 
ionization  exists  under  ambient  night  conditions. 

The  same  would  appear  to  be  true  at  first  sight  for  the  X  wave 
under  ambient  day  conditions  (Fig.  8) .  However,  this  is  somewhat  de¬ 
ceiving.  For  vertical  incidence  the  conduction  band  shown  by  vertical 
striping  in  Fig.  8  is  practically  transparent,  so  that  no  major  reflec¬ 
tion  of  the  X  wave  occurs  for  vertical  incidence  below  the  electronic 
gyrofrequency.  Nevertheless,  the  total  attenuation  of  the  transmitted 
X  wave  is  significant  except  well  down  in  the  hydromagnetic  band,  and 
it  becomes  prohibitive  above  about  104  Hz.  Thus  it  is  more  appropriate 
in  Table  1  to  list  the  inferior  penetration  frequency  for  the  X  wave 
at  vertical  incidence  under  ambient  day  conditions  as  about  104  Hz. 

Moreover,  for  the  X  wave  under  ambient  day  conditions,  glancing 
incidence  differs  from  vertical  incidence  in  a  rather  significant  way. 
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For  glancing  incidence  (i  =  80  deg)  the  curve  for  1  neper  attenuation 

4 

in  Fig.  8  intersects  the  top  of  the  conduction  band  at  10  -Hz  and  again 
at  2.5  x  10^  Hz.  Between  these  two  frequencies  reflection  of  the  X 
wave  at  glancing  incidence  would  therefore  take  place  from  the  conduc¬ 
tion  band  under  ambient  day  conditions.  It  should  be  noted  from  the 
right-hand  scale,  however,  that  attenuation  below  the  level  of  reflec¬ 
tion  is  prohibitively  high  above  about  10^  Hz. 

For  the  0  wave  under  ambient  day  conditions  (Fig.  7)  the  inferior 
penetration  frequency  of  the  ionosphere  occurs  where  the  wavelength  is 
long  enough  to  make  the  combined  conduction  and  stop  bands  less  than  a 
"skin  depth"  in  thickness.  This  occurs  at  about  1.6  Hz  for  all  angles 
of  incidence.  Even  major  disturbance  of  the  D  region  of  the  ionosphere 
(Fig.  9)  only  drops  this  to  0.4  Hz. 

At  night  major  disturbance  of  the  D  region  drops  the  inferior 
penetration  frequency  of  the  ionosphere  for  the  0  wave  from  roughly 
the  ionic  gyrofrequency  to  0.6  Hz,  a  value  intermediate  between  1.6 
Hz  for  ambient  day  conditions  and  0.4  Hz  for  disturbed  day  conditions 
(Table  1).  Thus,  for  disturbed  night  conditions,  as  well  as  for  day 
conditions  (ambient  or  disturbed) ,  the  inferior  penetration  frequency 
of  the  ionosphere  for  the  0  wave  is  in  the  general  vicinity  of  1  Hz. 

It  is  only  under  undisturbed  night  conditions  that  it  rises  towards 
the  ionic  gyrofrequency  for  molecules  (around  30  Hz) . 

The  inferior  penetration  frequency  of  the  ionosphere  for  the  X 
wave  is  much  more  dependent  upon  the  degree  of  disturbance.  Whereas 
for  ambient  night  conditions  the  inferior  penetration  frequency  of  the 
ionosphere  for  the  X  wave  is  not  much  less  than  the  electronic  gyro¬ 
frequency,  it  drops  to  a  few  hundred  hertz  under  disturbed  night  con¬ 
ditions.  It  drops  to  a  few  tens  of  hertz  under  disturbed  day  conditions 

4 

from  an  ambient  day  value  of  about  10  as  discussed  above. 

Curves  showing  the  overall  attenuation  for  transmission  through 
the  ionosphere  to  a  high  satellite  using  the  ionospheric  models  under 
discussion  are  given  elsewhere. ^  The  same  report  presents  curves 
for  the  variation  of  the  rate  of  attenuation  with  height  for  numerous 
frequencies  from  1  Hz  to  100  kHz.  A  feature  of  these  curves  at  fre¬ 
quencies  below  about  100  Hz  is  that  they  show,  in  addition  to  the  usual 
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electronic  behavior,  maxima  of  attenuation  in  the  vicinity  of  the 
level  where  (between  110  and  120  km)  arising  from  ionic  col¬ 

lisions  with  neutral  molecules.  These  maxima  are  of  particular  inter¬ 
est  in  the  neighborhood  of  the  ionic  gyrofrequency .  This  is  the  cause 
of  the  dip  in  the  curve  for  1  neper  penetration  near  the  ionic  gyro- 
frequency  in  Fig.  8  and  for  the  kink  in  the  same  curve  in  Fig.  10. 
Another  interesting  feature  of  the  results  in  Ref.  1  is  that  under  am¬ 
bient  conditions  there  is  actually  a  local  maximum  in  the  transmission 
loss  of  the  X  wave  through  the  ionosphere  near  the  ionic  gyrofrequency 
of  the  order  of  2  dB  by  night  and  11  dB  by  day,  virtually  independent 
of  angle  of  incidence.  Under  disturbed  conditions,  however,  this  local 
maximum  becomes  an  almost  insignificant  ledge  on  a  rapid  rise  of  trans¬ 
mission  loss  with  frequency  near  the  ionic  gyrofrequency. 
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APPENDIX 


In  this  Appendix  we  list  the  parameters  and  equations  used  to  ob¬ 
tain  the  numerical  results  given  in  this  report  and  in  Ref.  1.  The 
equations  are  the  well  known  ones  of  the  magneto-ionic  theory  includ¬ 
ing  positive  and  negative  ions  and  are  stated  without  derivation.  A 
discussion  of  the  range  of  applicability  of  the  equations  is  given  in 
Section  II.  We  begin  by  listing  parameters  and  defining  symbols  and 
then  give  the  equations: 

-3 

N  =  electron  density,  m 

e  -3 

N.  =  density  of  negatively  charged  atoms,  m 

la_  -3 

=  density  of  negatively  charged  molecules,  m 

N^a+  *=  density  of  positively  charged  atoms,  m  3 

t  =  density  of  positively  charged  molecules,  m  3 

0  =  temperature,  °K 

The  various  charged-particle  density  profiles  were  used  to  obtain 

Fig.  1.  The  ions  were  assumed  molecular  below  150  km  and  atomic  above 

160  km;  ©  was  obtained  from  Ref.  (13);  Z  »  p/1.225,  where  p  is  neutral 

-3 

atmospheric  density  in  kg  m  and  was  taken  from  Ref.  3.  The  electron- 
neutral  and  ion-neutral  collision  frequencies  vg  and  were  obtained 
from  the  simple  approximate  relations 


v  =  6.1  x  10  9  (^)~3/2  N.  +  1.8  x  1011  Z  sec 
e  V300/  i 

9  -1 

=  4  x  10  Z  sec 

Constants  which  were  used  are 

a  =  6370  km 

e  =  1.602  x  10-19  C 

m  *=  9.106  x  10~31  kg 

e  ° 

m  =  16  x  1837  x  m 
a  e 

m  =  29  x  1837  x  m 
m  o  e 

c  =  2.998  x  10  m/sec 

e  «  8.854  x  io"12  F/m 

o 
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The  squares  of  the  various  angular  plasma  frequencies  (in  hertz) 
are  given  by 


2  -  e  » 

WNe  e  m  e 
o  e 


2  -  e  XT 

WNa-  e  m  ia- 
o  a 


2  e  „ 

UNm-  t  m  im- 
o  m 


2  e  XT 

<!>„  _1_  =  -  N- 

Na+  e  m  ia+ 
o  a 


2  e  XT 

“XT  ,  “  -  N.  . 

Nm-t-  e  m  inH- 
o  m 


and  the  angular  gyrofrequencies  by 


or,  =  2tt  x  1.7  x  106(1  -  3-) 
Me  a 


Ma  16  x  1837  Me 


°Mm  29  x  1837  ^e 


OJ  =  2  TTf 


where  Z  is  the  height  above  the  ground  in  kilometers  and  f  is  the  wave 
frequency  in  hertz. 

We  now  list  the  following  functions,  where  the  upper  sign  corres¬ 
ponds  to  the  0  wave  and  the  lower  sign  to  the  X  wave. 
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A  = 


“Ne  (“  1  “Me5 
<*>[(«*>  ±  “Me)2  +  Ve] 


A 

a- 


“Na-  (“  1  t°Ma) 
“K“  ±  “Ma)2  +  vil 


m- 


“Nm-  (“  *  “Mm5 

“[<“  *  “W2  +  v±] 


^  V+  (“  T  t°Ma) 

a+  u[(o>  =F  O)^)2  +  V2] 


.  Vf  (ai  T  “Mm5 

tn+  ,  ,  v  2  ,  2. 

“[(“*“Mm5  +  Vi] 


B  = 


2 

“vr  V 

Ne  e 


“[(oJ  *  “Me)2  +  ve] 


2 

u)„  v . 
Na-  l 


a-  r/  .  \2  .  2, 

“[(“  ±  “(fe)  +  vi^ 


2 

Cl3.  _  V. 

Nm-  i 


m-  .  n2  ,  2, 

“[(“  4  “Mm5  +  Vi] 


<D„  ,  V  . 
Na+  l 


a+  r .  _  \  2  2 . 

T  “Ma5  Vi] 


to  v 
Nmf  i 


m+  r,  _  v  2  2, 

“[(“T“Mm5  +  vi] 
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A  =  A  +  A  +  A  +A.+A. 
e  a-  m-  a+  m+ 


B  =  B+  B  +  B  +  B  ^  +  B  , 
e  a-  m-  a+  m+ 


The  refractive  index  is  given  by 


n  =  (1  -  A  -  JB)- 


R(n)  >  0,  I(n)  <  0 


where  j  =  /-I ,  whence  the  local  wavelength  can  be  computed  from 


‘-fife-*10'3 


By  defining 


2  1/2 

q  =  (cos  i  -  A  -  jB)  R(q)  >  0,  I(q)  <  0 


a„,  =  8.686  —  I(q)  x  10^  dB/km 
T  c 


the  transmission  loss  is  computed  from 


T  =  j  aTdz  dB 


After  defining 


=  8.686  ^  B  sec  i  x  10^  dB/km 


the  reflection  loss  was  computed  from 


R  =  1  oi  dz  dB 

Jq  R 
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